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Abstract 

A self-consistent solution for a thin accretion disk with turbulent convection is presented. 
The turbulent convection plays a double role: it provides the disk viscosity and takes part 
in the vertical transport of the released energy. Rather than assuming arbitrary phenomeno- 
logical parameterizations for the disk viscosity, the latter is derived from a physical model 
for turbulence. Employing this model, we express the turbulent viscosity and the vertically 
averaged convective flux in terms of the local physical conditions of the disk which, in turn, 
are controlled by the former two. The resulting self-consistent disk structure, and the ratio 
between the convective and total fluxes are obtained for radiation and gas pressure dominated 
regions, and for electron scattering and free-free absorption opacities. 

In the gas pressure region, two distinct solutions are obtained. In one, the convective 
flux is much larger than the radiative flux and the blackbody region extends over the entire 
gas pressure region and could also extend down to the inner boundary of the disk. In this 
solution the temperature proflle is close to adiabatic. In the other solution, the convective 
flux is about a third of the total flux, the dimensionless superadiabatic temperature gradient 
is ~ 0.6 and there exist the gas pressure blackbody and electron scattering regions as well as 
the radiation pressure region. 

In the radiation pressure region, the temperature proflle is very close to adiabatic, and 
the disk is geometrically thin and optically thick even for super Eddington accretion rates. 
The fraction of the convective flux, out of the total flux, increases with the accretion rate, and 
for accretion rates comparable to the Eddington limit is close to 1. This variation stabilizes 
the, radiation pressure region, so that unlike the a disk, all the disk regions are secularily 
stable. The values of the effective a-parameter are rather small: < 5 X 10""^, ~ 1 X 10~^ and 
~ 5 X 10~^ for radiation pressure region and for the two solutions in the gas pressure region, 
respectively. 

Subject headings: accretion, accretion disks — convection — turbulence 
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1. Introduction 



There is quite a large body of observational data suggesting the existence of accretion 
disks around protostars, compact stellar objects and active galactic nuclei (AGNs). No won- 
der that accretion disk models are extensively employed to interpret these observations. How- 
ever, to a large extent the models used are more a descriptive than predictive tool (Pringle 
1981). The reason for this state of affairs is the lack of a physical model for the disk turbulent 
viscosity, i^^, usually parametrized in the form (Shakura & Sunyaev 1973) 

vt = acsh. (1) 

Here Cg is the speed of sound, h is a characteristic scale height, and a is a dimensionless 
parameter which neither value nor dependence on the physical conditions in the disk are 
known. Thus, often, the observations are used to fit the a-parameter for a given system. 
Usually a constant a is assumed for gas pressure dominated regions of the disk. For radiation 
pressure dominated regions both the former prescription and a proportional to the ratio 
between the gas and total pressure, have been suggested. In the lack of a physical model for 
the turbulent viscosity it is not clear which prescription, if any, is the correct one. 

A self-consistent solution to the disk equations requires that the turbulent viscosity be 
determined by the physical conditions of the disk. In turn, the resulting turbulent viscosity 
will control the same physical conditions. The determination of the turbulent viscosity re- 
quires a model for turbulence and the specification of the particular instability that generates 
the turbulence. The predictions of the self-consistent solution, obtained for the assumed gen- 
eration mechanism, can then be compared to the observations thus allowing to decide how 
relevant is the above generation mechanism to the astrophysical system under study. 

In spite the very large Reynolds numbers, Keplerian disks are stable against shear- 
generated turbulence, at least in the linear analysis (Pringle 1981). On the other hand it 
has been pointed out that they are unstable against turbulent convection, both in radia- 
tion pressure and gas pressure regions and for various opacities (e.g., Bisnovatyi-Kogan & 
Billinikov 1977; Shakura, Sunyaev, & Zilitinkevich 1978; Lin & Paploizou 1980). Therefore, 
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in this paper we implement the above proposed selfconsistent scheme for the case of turbu- 
lent convection in a thin accretion disk. The turbulence model employed (Canuto, Goldman 
& Chasnov 1987) provides vertically averaged values for the turbulent viscosity and for the 
convective flux. Thus, no attempt to resolve the vertical structure of the disk is made, and 
all quantities are either vertical averages or midplane values. A detailed vertical structure 
requires a turbulence model, that could yield the z-dependent values of the turbulent viscosity 
and of the convective flux. Such a model would be considerably more complex as the effects 
of the turbulence could no longer be represented by a turbulent viscosity and a convective 
flux but include also diffusion terms of various turbulent ensemble averages, which vanish in 
a vertically averaged model. In this work, we wish to focus on the more general features of 
the proposed self-consistent approach and avoid the complications of a z-dependent model. 
We note that turbulence bulk properties such as the turbulent viscosity and the convective 
flux are contributed mainly by the largest eddies. The vertical extent of the latter is ~ /i 
and thus vertical averaging is expected to represent fairly well their contribution. Indeed 
such averaging was shown to yield quite good predictions regarding laboratory convection 
(Canuto et al. 1987). Thus, our working assumption is that the important characteristics of 
the self-consistent solution would be evident also in the simpler vertically averaged approach. 

Doubts concerning the direction of the angular momentum flux associated with convec- 
tion have been raised by Ryu & Goodman (1992) and by Kley, Papaloizou, & lin (1993). In 
the flrst work linear convectives modes were shown to give rise to an inward angular momen- 
tum flux. However, as noted by the authors, this result may just reflect the inconsistency 
of their model and in a fully developed turbulence the situation would probably be differ- 
ent. Kley eA al. (1993) assumed an imposed underlying large viscosity and performed direct 
numerical simulation for convective perturbations. They interpret their results as indicating 
an inward angular momentum flux advected by convection. As the authors note this may be 
the result of the large value of the non convective viscosity that was used. Indeed, Cabot & 
Pollack (1992) flnd that when the Reynold numbers are increased ( viscosity lowered) tur- 
bulent convection provides a positive turbulent viscosity, thus ensuring an outward angular 
momentum transport. Since in the present model the turbulent viscosity is positive deflnite 
(see eqs. [13] and [19]) there will be an outward flux of angular momentum as required. 
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We have in mind disks surrounding compact objects of sizes ranging from stellar to 
galactic. For AGN accretion disks, the radiation pressure region can encompass a large part 
of the disk region from which most of the radiation is emitted. Therefore, we are interested in a 
self-consistent disk description for the following three regions: radiation pressure dominated 
region, gas pressure dominated region with electron scattering opacity, and gas pressure 
dominated region with free-free absorption opacity. Since turbulent convection is taken to be 
the source of the disk turbulent viscosity, the convective flux could play an important role in 
the transfer of energy to the disk surface. 

The resulting disk structure, and the ratio between the convective and total fluxes at 
given distance from the central object, are obtained as functions of S- the vertically averaged 
local value of the dimensionless superadiabatic temperature gradient. Practically, it is more 
convenient to express S in terms of i^o^ the ratio between the convective and total fluxes 
and use the latter in the parameterization of the disk solutions. Unlike the a-parameter, S 
could be determined once the vertical structure of the disk is obtained. Even in the vertically 
averaged approximation, considered in the present work, we derive general upper bounds on S. 
Moreover, using approximate vertically averaged relations we estimate the ratio between the 
surface and midplane temperatures. This together with an expression for the flux emanating 
from the disk, in terms of the surface temperature, provides an additional relation that can be 
used to flnd i^o and S. Due to the vertical averaging, the obtained solutions are correct up to 
factors of order unity. Nevertheless, we expect the results to be qualitatively representative 
of those of a detailed 2;-dependent model. 

Our main results are (1) In the radiation pressure dominated region the dimensionless 
superadiabatic temperature gradient is very small, so that the vertical temperature proflle 
is very close to adiabatic. This result is independent of the ratio of the convective flux to 
total flux. The ratio of convective to total flux is larger the larger is the accretion rate. For 
accretion rates comparable to the Eddington limit it varies between ~ 0.5 and close to 1. 
(2) For the gas pressure dominated region we flnd two distinct types of solutions. In one, 
practically all the flux is transported by convection and the dimensionless superadiabatic 
temperature gradient is small ( < 0.1). The blackbody region encompasses the entire gas 
pressure region and could extend down to a small distance from the central object, without 
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any radiation pressure region. In the other solution, the convective flux comprises about 1/3 
of the total flux and the dimensionless superadiabatic gradient is < 0.6. In this case all 
three regions of the disk exist. 3) The effective a-parameter is expressed as a function of the 
local dimensionless superadiabatic temperature gradient. Its typical values are quite small 
(this is a general feature of any three dimensional turbulence which has no energy sources 
external to the disk; Goldman 1991). 4) In the radiation pressure region an increase in the 
luminosity leads to an increase in the fraction of the convective flux which in turn leads to a 
net increase in the surface density. Consequently, the surface density is an increasing function 
of the luminosity, in all disk regions. Therefore, all regions are secularily stable, in contrast 
to a disks that are unstable in the radiation pressure region. 

The role of convection in accretion disks has been discussed by many authors within the 
phenomenological a -parameter model for the turbulent viscosity, and the phenomenological 
mixing length approach for the convective flux (see, e.g., Bisnovatyi-Kogan & Billinikov 1977; 
Livio & Shaviv 1977; Liang 1977; Shakura eA al. 1978; Vila 1978; Lin & Paploizou 1980; Tayler 
1980; Smak 1982; Ruden & Lin 1986; Duschl 1989, Milsom, Chen, & Taam 1994). 

Cabot et al. (1987a, b) applied an earlier, less developed, version of the turbulence model 
employed here (Canuto & Goldman 1985) to the protosolar accretion disk. The dominant 
opacity for that cool disk is due to dust grains. As stated above, we are interested here in 
much hotter disks surrounding compact objects, with radiation pressure dominated regions 
and with different opacities. Recently, Cabot et al. (1990) and Cabot & Pollack (1992) 
applied direct numerical simulations to convective disks with uniform and differential rotation, 
respectively. However, as is the case with direct numerical simulation in other problems 
involving turbulence, they are limited to low Reynolds numbers which are orders of magnitude 
smaller than those characterizing astrophysical accretion disks. 

Some of the ideas discussed here were presented in an earlier work (Wandel & Goldman 
1991). 

2. Disk Equations 

Let us consider a stationary geometrically thin accretion disk in which h <^ where 
h is the disk half- width at a distance R from the central object. The radial gradient of the 
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pressure is small compared to the gravitational force per unit volume exerted by the central 
object, so that the disk angular velocity is Keplerian 

O=(f^)'^'=2xl0-X-V-^/^3-\ (2) 

with M the central object mass, Mg = MI{WMq ), and r = Rc^/{GM). 

The vertically averaged disk equations include the energy equation, the angular mo- 
mentum equation, and the hydrostatic equilibrium equation. The energy released by the 
turbulent viscosity, per unit area of each face of the disk, Q, is given by (Pringle 1981) 

Q = ^iytphQ\ (3) 

where p is the density at the disk midplane and is the vertically averaged value of the 
turbulent viscosity. The angular momentum equation expressing outward radial transfer of 
angular momentum, due to the interaction of the turbulent viscosity with the Keplerian shear, 
is 

OM(^(r) = GTTflhput. (4) 

where M is the rate of mass accretion, and (f>{r) accounts for the inner, stress-free, boundary 
of the disk; for a nonrotating black hole (f>{r) = 1 — Combining equations. (3) and 

(4) results in 

Q = —fl^M^ir) = 1.2 X 10^°L^M-\--\Tgs-\ (5) 
Svr 

where equation (2) was used to express O, and where the accretion rate M was expressed in 
terms of the modified Eddington ratio , 

£. = A^,,,=0.4(^^)M-«r). (61 

Here Le is the Eddington luminosity, where the efficiency for release of gravitational binding 
energy was taken to be that of a nonrotating black hole. Equation (5) reflects the fact that 
the ultimate energy source is the gravitational binding energy, which is unlocked by the 
interaction of the turbulent viscosity with the Keplerian shear. In a stationary disk 
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F = Q. (7) 

with F denoting the energy flux emerging from each side of the disk surface. Since in this work 
we do not solve for a detailed vertical structure, the flux F is taken to be a vertical average. 
Since also Q, in equation (3), is represented in terms of vertically averaged quantities, we 
interpret now equation (7) to represent a relation between vertical averages. This is true up 
to factors of order unity that depend on the details of the vertical structure. 

The vertically averaged equation of hydrostatic equilibrium (up to a factor of order unity 
which depends on the detailed vertical structure) is 

P = pQ^h^, (8) 

where P is the total pressure at the disk midplane. As seen below, turbulent convection is 
subsonic so the turbulent pressure can be neglected compared to the thermal pressure, so 
that 



with the gas pressure, Pg 



and the radiation pressure, Pr 



P = Pg+ Pr, (9) 



Pg ^ (10) 



Pr = laT\ (11) 

where is the Boltzman constant, m is the atomic mass per particle (m = 0.59mp for 
solar abundances), a is the blackbody constant, and T is the midplane temperature. For 
convenience we introduce the parameter /3 — the ratio between the gas pressure and the 
total pressure 



(12) 



3. Turbulent Convective Viscosity 

The model for turbulence employed (Canuto et al. 1987) yields the turbulence spectrum 
in terms of n{k) — the rate controlling the energy input, at wavenumber A;, from the gen- 
erating source into the turbulence . The turbulent viscosity, which is contributed by all the 
eddies comprising the turbulence spectrum, can nevertheless be expressed simply as 

vt = n{ko)k^'^, (13) 

where ko is the smallest wavenumber (largest eddy) present in the turbulence spectrum. In 
principle, n(k) depends also on the turbulence spectrum itself. However, Canuto et al. (1987) 
demonstrated that in the case of turbulent convection, it can be approximated by the growth 
rate of the unstable modes of the linearized equations. These authors found that the resulting 
turbulence spectrum and convective flux are in good agreement with those derived from a 
self-consistent renormalized n(A;), and with laboratory experiments on convection at high 
Rayleigh numbers. 

Turbulent convection in accretion disks is likely to be affected by the interaction of the 
strong Keplerian shear with the eddies. In the turbulence model, employed here, the effects 
of the shear are introduced through the growth rate n(k) which corresponds to convection 
in a sheared (differentially rotating) Keplerian disk (Goldreich & Schubert 1967; Canuto, 
Goldman, & Hubickyj 1984). The growth rate for the largest eddies (smallest wavenumbers) 
which determines the turbulent viscosity in equation (13), can be expressed in terms of a 
dimensionless growth rate N multiplied by the inverse of the buoyancy timescale 

/ 2^x1/2 

n(ko) = N ig.aS-j = NQS^^^ (14) 

where N is the positive real part of the solution of a third order algebraic equation (Canuto 
et al. 1984). Its dependence on the physical parameters will be considered in § 5. In equation 
(14), gz = Q'^ z is the vertical gravitational acceleration at height z due to the central object, 
a is the coefficient of thermal volume expansion at constant pressure (Chandrasekhar 1967) 
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_ 1 fdV\ , (4-3/3) 

The dimensionless superadiabatic temperature gradient S is given by 

-(-H(f)-(f)J)' 

with the angular brackets denoting vertical average. The parameter S appearing in the right 
handside of equation (14) is given by 

* = ^i^*. (17) 

For gas pressure dominated regions 8 = 8 while for radiation pressure dominated regions 8 = 
4:fi~^8 ^ 8. All the quantities in equation (14) represent vertical averages. The wavenumber 
is expressible in terms of the disk half- width h 

ko = h-^Ti{l + xfl^, (18) 

where x = k'^/k'j, is a measure of the anisotropy of the largest eddies, kp and k^ are the 
horizontal and vertical wavenumbers respectively, so that k'^ = k^ -\- k^ with k^h = vr. 
Using equations (14) and (18) in equation (13) yields 

Before discussing the dependence of N on the physical parameters, we present the expression 
for the convective flux. 

4. Turbulent Convective Flux 

Since in the present case convection is the generating mechanism of the disk turbulence, 
it is conceivable that the convective flux could be important in the transfer of the released 
energy to the disk surface. Therefore, in equation (7) the vertically averaged flux is taken to 
be the sum of the radiative and convective fluxes 

F = Fr + F, . (20) 
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The radiative flux Fr is given by 



4acT4 . X 

F. = — , (21) 

with T = uph the optical depth (T, /?, and the opacity n are midplane values). The vertically 
averaged convective flux, Fc is given by (Canuto eA al. 1987) 

Fc = Cpp^A^^n\U)k-^, (22) 

with A a dimensionless constant that depends on the strength of convection, as detailed in § 
5. The speciflc heat at constant pressure Cp is (Chandrasekhar 1967) 

r2 A;b (4-3/3) P V2 _^2,2 ...^ 
Cn = r:^ = a — = — ail n , zo 

^ r2 - 1 m /32 r2 - 1 p r2 - 1 ' ^ ^ 

where denotes the second adiabatic index — equation (8) has been used to obtain the 
third equality. From equations (14), (15), (18), (22), and (23) follows that 

F, = — ^-^ph^9.^NW\ (24) 

The fraction of the convective flux out of the total flux, 1^0, can be obtained from equations 
(3), (19), and (24): 

Before applying equations (19) and (24) to the disk equations, we consider below the depen- 
dence of the dimensionless growth rate, iV, on the physical parameters of the disk. 

5. The Limits of Strong and Moderate Convection 

Using the standard expressions for the kinematic radiative viscosity (Weinberg 1971) and 
for the kinematic plasma viscosity (Spitzer 1962) we flnd that, for the physical conditions 
typical to accretion disks, both are much smaller than the radiative conductivity. Thus, the 
Prandtl number (the dimensionless ratio between the kinematic viscosity and the radiative 
conductivity) is very small compared to unity. Only for an extremely high radiation pressure 
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is the radiative viscosity comparable to the radiative conductivity. In the hmit of very smaU 
Prandtl number, the dimensionless growth rate N depends only on the value of S — the 
dimensionless product of the Rayleigh and the Prandtl numbers (Canuto et al. 1984). For a 
thin disk 

S = g^aSjh\-^ = n^hHx~\ (26) 
where x is the radiative conductivity 

F h 

X = (27) 

The possible solutions of the equation for N{S) (Canuto eA al. 1984) can be classified according 
to the strength of the convection, measured by the dimensionless parameter S"*, 

^ - vr2(l + x)- ^^^^ 
Strong convection is obtained in the limit S"* ^ 1, weak convection corresponds to S"* ^ 1 
and moderate convection to S* ~ 1. Using equations (19) and (26) in equation (28) one finds 

S* = -, (29) 

X 

which has a simple physical interpretation as the ratio between the timescales characterizing 
the radiative transport and the turbulence. Using now equation (3) to express Vt in terms of 
Q and equation (27) to express x ii^ terms of i^^, as well as equations (15) and (23), one finds 

4 (4-3/3) Q 4 (4-3/3) , 

^ =9r.-i /3 ^ = 9r.-i /3 ^^-^"^ ' 

where i^o is the ratio of convective to total flux, see equation (25). Equation (30) must 
be satisfled in any self-consistent disk solution in which the turbulent viscosity is due to 
convection. In radiation pressure dominated regions (/3 ^ 1) strong convection is guaranteed, 
even if the convective flux were small compared to the radiative flux (small i^o)- In gas pressure 
dominated regions(/3 — )• 1) equation (30) yields S* = 1.1(1 — i^o)~"^ , for V2 =5/3. Thus, in this 
case convection is either moderate or strong depending on the value of 1^0 • Strong convection 
occurs only for 1^0 close to unity, namely only when the convective flux is much larger than 
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the radiative flux. In each of these two hmits we would express N and x in terms of (5, so 
that both the turbulent viscosity and the convective flux will depend only on S. Practically, 
it will be more convenient to use equation (25) to express S in terms of i^o, and use the latter 
as the parameter in the various expressions. 

5.1. Strong Convection 

In the strong convection limit one has (Canuto et al. 1984) 



For a given superadiabatic gradient there are many possible unstable modes, correspond- 
ing to different values of the anisotropy parameter x, each resulting in a different value for 
the turbulent viscosity. Generally, the latter will be a superposition over the different x val- 
ues. Since the disk structure will be dominated by those modes that contribute most to the 
turbulent viscosity and hence to the energy production, we approximate the superposition 
by that x which maximizes the turbulent viscosity for a given superadiabatic gradient. From 
equations (19) and (31) then follows that 



x = 1(1 +3^1) (32) 



for which equation (31) yields 



N = 3"^/^ (33) 
so that the turbulent viscosity, equation (19), becomes 

z/i = 0.039 =0/i^ (34) 

In the strong-convection limit we flnd, using equation (22) with equation (52) of Canuto 
eA al. (1987), that the parameter A appearing in the expression for the convective flux, 
equation (24) is ~ 3. For the above values of N and A, equation (25) yields 



(35) 



which for radiation pressure dominated regions (F = 4/3) is 



Co = 17.6^, (36) 
while for gas pressure dominated regions (F = 5/3), 



Co = n.os. (37) 

Since 1^0 ^ 1 by definition, it foUows that for radiation pressure S < 0.057 and for gas pressure 
S < 0.09. Thus, the temperature gradient is close to adiabatic; in particular so in the case of 
radiation pressure for which S <^ S. Such a situation is indeed expected for strong convection 
which acts to reduce the temperature gradient. The corresponding values of x are quite large: 
X > 25 for radiation pressure and x > 17 for gas pressure dominated regions. This implies that 
the largest turbulent eddies are anisotropic with the horizontal dimension smaller than the 
vertical dimension by a factor of x^^'^ ( > 4 — 5). Such anisotropy, resulting from the Coriolis 
force, is a general feature of any three dimensional turbulence in a differentially rotating disk 
(Goldman 1991). 

Equation (34) can be used to define an effective a parameter 

'-■< 6''"- 

From equations (36) and (37) follows that cte// 5 5 X 10""^ for radiation pressure dominated 
regions and < 10~^ for gas pressure dominated regions. The small values of cte// result from 
the small value of the turbulent velocity compared to the sound velocity (see eq. [44])) as 
well as from the large anisotropy of the eddies with the horizontal dimension much smaller 
than the disk half-width. 



5.2. Moderate Convection 

As noted above, in radiation pressure regions convection is always strong but in gas 
pressure regions it can be also moderate. From Canuto et al. (1984) we find a solution for N 
for which S"* ~ 1 
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Therefore, the maximal turbulent viscosity is obtained for an anisotropy parameter 



x~^(l+2^-^) (40) 

yielding 

N ~ (18)"^/^ ~ 0.24. (41) 

Using equations (52) and (56) of Canuto et al. (1987), we find that for moderate convection, 
the dimensionless parameter A ~ 1. For the above values of A and N and for V = 5/3, 
equation (25) results in 

Co ~ 0.6(5. (42) 
Use of equations (40), (41) and (42) in equation (19) yields the turbulent viscosity 

vt^QMQ-^ =O/i^=0.05^^ — -Vth^. (43) 

1 + 8 1 + 2.5^0 ^ ' 

We note that moderate convection is obtained when S* is of order unity. Therefore, all the 
numerical values, in the relations above, should be regarded as representative, and could 
actually differ by factors of order unity. 

Equation (42) implies that 8 could be larger than in the case of strong convection (for 
gas pressure regions 8^8). For such larger values of 8, the anisotropy parameter, given by 
equation (40), is significantly smaller than in the case of strong convection. The effect of the 
smaller anisotropy parameter overcomes the decrease in N so that the effective a- parameter 
can be somewhat larger than in the case of strong convection. For example, i^o = 0.3 yields 
8 ~ 0.5, implying x ~ 2.5 and cte// ~ 5 X 10~^. 

6. Disk Structure — Strong Convection 

The expressions for the turbulent viscosity and for the turbulent convective flux, in the 
strong-convection limit, are used below to solve for the disk structure. This is done for 
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three different regions: radiation pressure dominated, gas pressure dominated with electron 
scattering opacity and gas pressure dominated, free-free absorption opacity, blackbody region. 
Given the mass of and the distance from the central object, and the accretion rate, we obtain 
the disk half-width /i, the midplane values of the temperature T and the density /?, and 
the ratio of the convective to total flux, as functions of the dimensionless superadiabatic 
temperature gradient. 

Before doing so we note that the turbulent pressure need not be included in the hy- 
drostatic equilibrium equation since the turbulent velocity is much smaller than the thermal 
velocity. The turbulent velocity (rms value), obtained using Canuto et al. (1987), is 

N 



vt ~ [A7r^n^{ko)k-^] 1/2 = A'/^ ——-—S'/^Qh, (44) 

(1 -|- x)'-'^ 

implying that the ratio of the turbulent to the thermal velocity (~ Qh) is < 0.06 for both 
radiation pressure and gas pressure dominated regions. 

6.1. Radiation Pressure Region 

In this region, /3 <^ 1, k = k^s = 0.4g~"^cm^, and the strong convection limit applies 
regardless of the value of i^o- Expressing the total pressure as 

P = -aT^^—-, (45) 
3 1-/3' ^ ^ 

the hydrostatic equilibrium equation, equation (8) takes the form 

^aT^^ = pQ'h\ (46) 

We use equation (25) to express the convective flux in terms of i^o and Q, so that equation 
(20) becomes, 

F, = (l-eo)Q, (47) 
which upon using equations (5) and (21) yields 

= (1 - i,) — ^^M<i>{r) = 1.2 X 10^^(1 - i,)L.M-\-\rgs-\ (48) 
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The angular momentum equation, equation (4) with Vt given by equation (34) resuhs in 

^3/2 .3/2 

MSir) = 0.753pO/i^ = = O.OlpO/i^ ^ -, (49) 

^ ' 1 + 8 ' 1 + 0.057^0 

where the second equahty foUows from equation (36). In what foUows the denominator in 
the right hand side of equation (49) wiU be approximated by 1 (since i^o ^ !)• 

Equations (46), (48) and (49), with equation (6), can be solved to obtain T, p and h 
as functions of i^o and the distance from the central object. Equation (36) can be used to 
express i^o in terms of 8. Combining equations (46) and (48) yields 

h = {lx W^cm) MsL, (1 - /3)"' (1 - Co), (50) 

so that 

| = 6.63i:,(l-/3)-^(l-eo)r-^ (51) 

Substituting h from equation (50) and O from equation (2), in equation (49) yields the 
midplane density p 

p = (1.2 X lO-^'gcm-') L;' (1 - (3f (1 - Cor'C^'^'r'/^ (52) 
which when used with equation (50) yields the optical depth for electron scattering 

Te, = K,,ph = (1 - (3f (1 - Cor'Co'^'r'/\ (53) 

The midplane temperature T is now obtained from equations (48) and (53) 

T = (3.7 X lO^'K) Mg"'/^ (1 - /3)'/' (1 - eo)"'/^eo"'^'^"'^' , (54) 
Equations (10), (11), (52), and (54) yield 

^ = 7.6 X 10^(1 - Co f/'fo^'LlM'J^ (1 - /3)-'/' r-''/\ (55) 

Therefore, a necessary condition for the radiation pressure region to exist is a value of 
1^0 not too close to unity or to zero. Thus, demanding /3 = 0.2 for r = 20 implies 2 X 
lO-iiL^^'/'Mg"'/' < Co < 1 -0.016i:;'/'Mg"'/'. 
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Formally, for i^o ^ 1 the above expressions for /i, /?, Tg^, and T, have the same form 
as in a a-disk with a small value of cte// (resulting from the small i^o)- Unlike in the a- 
disk model, the effective value of a is determined directly by the local physical conditions. 
However, as shown in § 8.5, in the radiation pressure region For i^o 1, the present 

convective disk is geometrically thinner and optically thicker the larger i^o is. Note that for 
1^0 — ^ 1, equation (51) implies that the disk is geometrically-thin at small r values even for 
super-Eddington accretion rates, L*^(r)~"^ > 1. An increase in i^o causes also an increase of 
T but the dependence is quite weak, see equation (54). In §8.5 we present an estimate for 
the dependence of i^o on the local disk parameters. 

6.2. Gas Pressure Region with Electron Scattering Opacity 

We consider first the limit of strong convection which in gas pressure dominated regions 
requires that i^o ^ 0.9, i.e. a convective flux which is at least an order of magnitude larger 
than the radiative flux. Now, n = k^s but 1-/3^1 and it is convenient to express the total 
pressure as 

P = ^pT, (56) 

mp 

which, when combined with equation (8), results in 

-T = Q^h^. (57) 



kB 



Equation (48) is unchanged in the present case. Using equation (37), the second equality in 
equation (49) becomes 



M6{r) = 0.02pO/i^ — ^ (58) 

^ ' 1 + 0.09^0 

and, here too, the denominator will be taken as equal 1. Combining equations (48), (57) and 
(58), with equation (6), yields the midplane temperature T, 

T = (1.5 X lO^JO M^^'^lVfi^lHl - eo)^/'eo"'^'^"'/'*'- (59) 
Further combination of equations (57) and (59) yields the half width of the disk 



3/2 



h = (7.3 X lO^'crn) Ml'''Ll'''r^l\l - (60) 

so that 

4 = 4.9 X IQ-^'M^^'^'lVfr^lHl - i. fl^'^i-^'^'r^l^^ (61) 
R 

The electron scattering optical depth is obtained from equations (48) and (59) 

Te. = 1.3 X IQ^VMl'^fi^'Hl - i, Vi^'''r-'l'. (62) 
The midplane density p is obtained directly from equations (60) and (62) 

p = (4.5 X IQ-^gcm-^) L'^'M"'/' V'/'(l " i,)-^'^' i-^^^'^'r-^^'^' . (63) 

The disk is geometrically thinner, optically thicker and has a lower midplane temperature 
the larger is i^o- 

6.3. Gas Pressure Blackbody Region 

In this region the dominant opacity is due to free-free absorption 

Kff = kopT-'l\ (64) 

where 

ko = 6.45 X 10^^ Cbfg~^cm^K^^^ (65) 

and Cbf is the bound-free enhancement factor which equals ~ 30 for solar abundances. 

As in the case of electron scattering opacity we consider first the strong-convection limit. 
Equations (56), (57) and (58) are unchanged in this case. The only change is in equation (48) 
where the opacity is now Hff instead of Hes- Equations (57) and (58) can be use to express 
h and p in terms of T, so that Hffoi equation (64) can also be expressed by T. Defining 

/r \ 1/20 

V~30" j ' ° 

equation (48), with Hff instead of now yields 
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T = (2.6 X lO'K) LV^''M;^'''^^'^il,r-^'\ (66) 
Substituting this T into equation (57) yields 

h = (3.0 X lO^'cm) Ml'^'Ll'^'r^'^ibbr'''' (67) 

so that 

4 = 2.0 X lQ-'M;''''Ll'''fr'l'i,y\ (68) 
R 

FinaUy, equations (58) and (67) resuh in 

p = (0.21,cm-^) M;'"'lV"'P'^I^^-:^-"\--''I^ (69) 

and 

Te. = 2.8 X lQ'Ml''LV''fi'l^i-,^i-''\-'l^ (70) 

7. Disk Structure — Moderate Convection 

As already noted, contrary to the radiation pressure region in which convection is always 
strong, in the gas pressure region convection can be either strong or moderate. The first 
possibility was addressed in § 6. In what follows we consider the second possible class 
of solutions for the gas pressure region — that of moderate convection. In this case too, 
equation (44) implies that the turbulent velocity is much smaller than the thermal velocity 
thus the turbulent pressure need not be taken into account. 

Equations (48) and (57) are unchanged while Vt is now given by equation (43) and i^o 
is related to 8 by equation (42). Therefore, instead of equation (58) the angular momentum 
equation will be 
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where equation (42) was used to obtain the second equahty. We wish to stress again that 
since for moderate convection S"* ~ 1 up to a factor of order unity, the relation between i^o 
and S (eq. [42]) as weU as equation (71) can also vary by factors of order unity. 



7.1. Gas Pressure Electron Scattering Region 

Repeating the same steps as in §6.2., and defining 



e* = (i-eo)'/'eo"'^'''(i + 2.5^0)^/' (72) 



we obtain now 



T = (6.9 X 10' K) Mg"'/'i:'/'/3i/^e*r-'/'°, (73) 



h = (4.9 X lO^^cm) M!^''L'J'r'^'U'^'r^'^'\ (74) 



4 = 3.3 X 10-3M-'/'*'i:y'/3-2/5e,i/Vi/2o, (75) 
R 



Te. = 5.7 X 10' lTm'J' (3'/' (1 - eo)-'d^-^/^ (76) 

p = (2.9 X lO-^'gcm-') lV M;' ' (i' ' {1 - i^yHl'^r-^^'^' ■ (77) 

7.2. Gas Pressure Blackbody Region 

Repeating the same steps as in § 6.3. we obtain now 

1/10 



T = (1.2 X 10^ JO [i^l^ ( ^ j L^/^^M-^/'e 



.r-'l\ (78) 



/n \ 1/20 

h = (2.0 X lO^'cm) Ml"'Ll"'r"' (^) il"r- 



(79) 
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- = 1.35 X 10-'M-'/''Ll/''r'^' ^ ) eJ'r'/^ (80) 



p = (1.4 X lO-'gcm-') M;" lV fi'^ I ' (1 - i,)-^lHV" r-^"\ (81) 



re. = 1.1 X lO^Mg^/^Ll/^V'/^ ^ ] ' i^-i^)-"'itr-"'. (82) 



-1/10 



8. Determination of i^o in the various Regions 

The solutions in the various regions of the disk depend on i^o — the fraction of the con- 
vective flux out of the total flux. This is determined in each region by 8 — the superadiabatic 
temperature gradient. Obviously, the latter two are known once the detailed vertical struc- 
ture is known. This however requires (see §1) a considerably more complex turbulence model 
than the one considered here. Thus, we wish to obtain an approximate determination of (5, 
and consequently of i^o within the framework of the vertically averaged approach. First we 
approximate the vertical average of the dimensionless temperature gradient by 

z dT\ hT — Tg Tg , . 

T^/ ~T^^ = ^- ^^^^ 

with Tg denoting the surface temperature at the disk top, z = h. Next, the vertically 

averaged dimensionless adiabatic temperature gradient can be approximated, making use of 

the hydrostatic equilibrium equation (8), as 

z dT\ V2-I 1 

T dz / V2 V2 

Thus, from equation (16) it follows that 
or equivalently 

T 1 

8 (86) 



- S 
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We note that, having been derived by an approximate verticaUy averaging, equations (85) 
and (86) are probably correct only up to factors of order unity. Even so, we expect the results 
to be representative of those corresponding to a more complex z-dependent modeling of the 
turbulence. 

Above the top of the convective disk, there is a surface layer (of width Ah << h) in which 
the flux is purely radiative. We further assume that this layer is isothermal at a temperature 
Tg. If this layer is optically thick to absorption, or absorption modifled by electron scattering, 
then the flux emerging from the disk can be represented by a blackbody or modifled blackbody 
with a temperature Tg. In this case equation (86), together with the disk equations, would 
enable the determination of all the disk variables, including S (or equivalently i^o)- In what 
follows, we implement the above procedure to the various regions of the disk. 

8.1. Gas Pressure Blackbody Region — Strong Convection 

In the strong convection limit S and i^o are related by equation (37). The latter relation 
together with equation (86) yields, for T2 = 5/3 

^-0.6-0.09^0 (87) 

In the blackbody region, Tg is determined by assuming a blackbody spectrum for the 
flux emerging from the disk surface 

T,= ( l^y^' = 1.2 X lO'K lI'^M;^'\-^I\ (88) 
\ac J 

where equation (5) has been used to obtain the second equality. Dividing this Tg by T of 
equation (66), and substituting in equation (87) yields an equation for 1^0 (which is equivalent 
to an equation for 8) 



(0.6 - 0.09^0) (1 - i^f'^'i^'''' = 4.6 X lO-^M-^/'^L,-^/' V"'/' 
resulting in 

Chf\ 



-1/10 

• (89) 



1 - Co ~ 3.8 X lO-i^M-^/^L.-^/^/J-^/^ (^^) , (90) 
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Thus, the convective flux is much larger than the radiative flux. With this value of i^o, we 
obtain from equation (37) that S = 0.09, and equation (87) yields Tg/T = 0.51. Substituting 
the above value of i^o in equations (66)-(70) results in 



T = (2.4 X lO^'K) lV^M;^'\-^I\ (91) 

4 = 6.1 X lO-^M-'/'Ly^-i/V^/^, (92) 
R 

T,, = 3xlO''M'J'LV'(3r-'/\ (93) 

, = (8.1^/cm-3) M-'/'L^/^^/V-i^/^. (94) 



This disk solution would have been obtained had we started with pure convective transfer 
(i^o = 1) and employed equations (56), (57), (58), (87), and (88). 

The line marked BBS in flgure 1 shows the L — t^s relation for the strong convection 
limit in the gas pressure, blackbody region. The positive slope implies that this solution is 
secularily stable. 

In order to flnd the boundary between the blackbody and electron scattering regions we 
evaluate the ratio between the free-free and electron scattering opacities. From equations 
(64), (65), (91), and (94) follows that 

^=2.1x 10' ( ^] mI''L-"'P'I\'I\ (95) 

Kes V 30 / 

which is larger than unity for all values of r. Therefore, it is possible for this solution to 
encompass the entire gas pressure region. As will be shown in § 8.2, there is indeed no gas 
pressure dominated electron scattering region, in the strong convection limit. 

In a-disks there is a radiation pressure dominated region for small values of r, when 
-L* ~ 1. Could the entire disk, in the present case, be a gas pressure dominated blackbody? 
To flnd out we compute the ratio between the gas and radiation pressures, in the gas pressure 
blackbody region. From equations (10), (11), (91) and (94) we obtain 
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^ = 3.25 X lO^L.-'/'Mg'/V'/'r^/'- (96) 

Pr 

Therefore, the disk could be gas pressure dominated blackbody for aU values of r, even when 
-L* ~ 1. Note that unlike in the a-disk, an increase of the mass of the central object favors 
gas pressure dominance. 

8.2. Gas Pressure Electron Scattering Region — Strong Convection 

In the electron scattering region, we apply the modified blackbody approximation to the 
emerging surface flux (Rybicki & Lightman 1979) 

which is valid when ^ ^^ff- The absorption opacity Hff is evaluated by substituting in 
equation (64) the surface temperature and the midplane density. One obtains from equations 
(5), (63), (64) and (97) the value of Tg which when divided by T of equation (59), yields 



^ = 2.8 X 10-^/3-^/1^ (^^^ i:;'/''M-'/''(l - eo)-'/''eo^''^"'/''- (98) 
Using this in equation (87) yields an equation similar to equation (89) 



(0.6 - 0.09^0) (1 - eo)'/''eo"'^'' = 2.8 x IQ-'r'^'' (f^) L^'/^'M^'/^'r-'/''. (99) 
Equation (99) yields a value of i^o which is extremely close to 1 

-17 A/r-2/3r -1/3/^-7/2 f£M.\ .-1/2 



1 -Co ~ 1 X 10-i^Mg-^/U,-^/^/3-^/^ (^^J r-i/^ (100) 

The ratio Hff/nes^ is obtained from equations (64) and (65), with T and p given by 
equations (59) and (63), respectively, and i^o given by equation (100) 

^ = 4.6 X 10^1 ( ^y^' M'J'L:'/'^^r\ (101) 
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The above ratio is larger than unity for aU values of r. This result is inconsistent with 
the definition of the electron scattering region. Since we are dealing with a thin disk, it 
is of interest to repeat the above procedure with Hff in equation (97) evaluated in terms 
of the mid-plan rather than the surface temperature. Doing so yields again i^o extremely 
close to unity, and the conclusion remains unchanged. Thus, we conclude that in the strong- 
convection limit there is no gas pressure dominated electron-scattering region and the entire 
gas pressure region is a blackbody, described by equations (91)-(94). 

8.3. Gas Pressure Blackbody Region — Moderate Convection 

We follow here the same procedure as in the case of strong convection, § 8.1. Equation 
(86) for T2 =5/3 yields, upon using equation (42) to express S in terms of 1^0, 

^-0.6-1.67^0 (102) 

From equations (72), (78), and (88) we obtain the ratio of T^/T, which when substituted in 
equation (102) results in 



(0.6-1.67eo)(l-eo)'/'''eo"'^'''(l + 2.5eo)^/' =0.1M-^/^''i:;^/^V-'/' [j^j . (103) 

For given values of the various dimensionless parameters, equation (103) can be solved 
for 1^0- Substitution of this 1^0 into equations (78)-(82) will yield the values of the disk physical 
variables. However, it is not possible to obtain an analytical expression for 1^0 in terms of the 
dimensionless parameters, as was done in § 8.1. There, it was possible because 1^0 was very 
close to unity, for all relevant values of the dimensionless parameters. 

The right hand side of equation (103) is insensitive to the values of the dimensionless 
parameters. For a representative case where they are set equal to unity the solution of 
equation (103) is 

Co ~ 0.32 (104) 
corresponding to (5 ~ 0.53 < l/r2. Substituting this 1^0 in equations (78) - (82) results in 
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T = (1.8 X 10^ JO ( ^ ] ' lI'^'m;^'\-^I\ (105) 

30 



1/10 



h = (2.5 X lO^^cm) Mj/^^L^y^V"'/' (f^) '^'^ ' (^O^) 



4 = 1.7 X 10-'M-'/''LV''r'^' ( ^V^'^r^/^ 
ii \ 30 y 

-3/20 



(107) 



p = (7.4 X lO-'gcm-') M;" lV I ' r-''/\ (108) 



Te, = 7.4xlO«M8'/'i::/'V/'(^^j ' (109) 

The line marked BBM in figure 1 shows the L — t^s relation for the moderate convection 
limit in the gas pressure, blackbody region. The positive slope of the curve implies that this 
solution is secularily stable. 

To find the boundary between the electron scattering and blackbody regions, as deter- 
mined in the blackbody region, we apply T and p of equations (105) and (108) to equation 
(64) and obtain 

^ = 1.4 X 10-2 ( ^] ^'^ L-"'p'l^r'l^ (110) 

Kes V 30 / 

implying that the blackbody region {nff> Hes) exists for 

r>2S^{^^y''\V'r"\ (111) 

but see equation (120). 

8.4. Gas Pressure Electron Scattering Region — Moderate Convection 

Similarly to the case of strong convection, the surface temperature obtained from the 
modified blackbody approximation, equation (97), and the temperature given by equation 
(73) are substituted in equation (102), resulting in 
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(0.6 - 1.67^0) (1 - eo)'/''eo"'^''(l + 2.5^0)''/'' = 

1.1 X 10-' r'^'' (f^) L^'/^'M^'/^'r-'/''. (112) 

Once the values of the various dimensionless parameters are specified, equation (112) 
can be solved for i^o- Substitution of this i^o into equations (73)-(77) will yield the values 
of the disk physical variables. However, it is not possible to obtain an analytical expression 
for in terms of the dimensionless parameters, as was done in §8.2. There, it was possible 
because i^o was very close to unity, for all relevant values of the dimensionless parameters. 

Setting all the dimensionless parameters equal to unity, equation (112) yields 

Co - 0.357 (113) 

corresponding to (5 ~ 0.59 < 1 /T2- We note that an increase of the right hand side of equation 
(112) by a factor of 10 yields i^o ~ 0.33, while decreasing the right hand side results in i^o 
larger than that of equation (113) but smaller than 0.36. For i^o = 0.357 equations (73) - 
(77) become 

T = (9.8 X 10' JO M-'^'LV'^'/'r-'/'\ (114) 

h = (5.8 X lO^^cm) Mj/'^Ly^/J-^/V^i/^", (115) 

4 = 3.9 X lO-^M-'/'^Ly'/J-^/V^/^", (116) 
R 

Te, = 3.6 X 10«i:'y'Mg'/'/3^/V-3/^ (117) 

p = (1.6 X 10-2^/cm-3) L'y^M-'/'V'/'r-^^/^". (118) 

The line marked ESM in figure 1 shows the L — t^s relation for the moderate convection 
limit in the gas pressure, electron scattering region. The positive slope implies that this 
solution is secularily stable. 
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To find the boundary between the electron scattering and blackbody regions, we substi- 
tute T and p from equations (114) and (118) to equation (64) and obtain 



""ff - 8.3 X lO-*^ L:^I3^I\^I^ ■ (119) 



implying that the electron scattering region {nff< Hes) exists for 



r < 2.4 X 10^ ■ LV^r^l\ (120) 

In determining the boundary between the electron scattering and the blackbody gas pressure 
regions both equations (HI) and (120) should be considered. 

In order to find the boundary between the electron scattering gas pressure region and 
the radiation pressure region we evaluate the ratio between the gas and radiation pressures in 
the gas pressure electron scattering region. Employing equations (10), (11), (114) and (118) 
results in 



-2/3 



^ = 9.4 X lO-^L.-'/'M-'/'V'/'r^'/'"- (121) 



Thus, Pg > Pr for 



r > imLl'''^^Ml'^^(i^l\ (122) 



8.5. Radiation Pressure Region 

In the radiation pressure dominated region equation (36) implies that 8 « 8 < 0.057, 
meaning that the disk is effectively adiabatic so that p . The density at the disk surface. 
Pa is expressible as 

Pa = p(Ta/Tf (123) 

In order for the modified blackbody approximation to apply, it is necessary that the 
effective optical depth for absorption of the surface layer 

/ \ 1/2 

r„_^, = /€e.p.A/i (124) 
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be larger than unity. Here Ah (<< h) is the width of the surface layer, and Hff is evaluated 
in terms of pa and Tg. 

In the surface layer the flux is purely radiative and constant. This means that the gra- 
dient of the radiative pressure divided by the density is constant. On the other hand the 
gravitational acceleration is proportional to 2;, thus only the (much smaller) gas pressure gra- 
dient divided by density is varying. Solving the set of coupled equations for the temperature 
and density we flnd that indeed the surface layer in this case is nearly isothermal, the density 
decreases as a gaussian, and Ah ~ h(l3Ta/Ty/'^ . Using this Ah and equations (52), (53), 
(54), (64), and (123), equation (124) yields 

r„,^. = i.ii:-/-M-^/-(i-/3)^/^*' (^y^"(i-eo)--/^''eo-"/'v-/- {f^^T- 

(125) 

Assuming for the moment that indeed this r^^.^, , > 1, we employ equation (97) with equations 
(52), (54), (64), and (123) to obtain 

^ = 5.14(1 - ^)-'/'' (^^^ lV''mI/'\i - i,f'l^'il''''r-'''''- (126) 
from which 

8^33 20^11 

(1 - eo)eo"/'' = QmM;"''L:''"'{l - fif^l'' r^'l^' (127) 

At each given radial distance, r, this constitutes an equation for 1^0 as function of the 
accretion rate (expressed via -L*), the mass of the compact object and the local physical 
conditions. The latter are manifested through /3, Chf and most importantly the ratio Tg/T. 
The estimate of equation (86), for T2 = 4/3, yields Tg/T ~ 0.75. As already noted, this is 
only an order of magnitude estimate. For this value, and taking r ~ 20, there is a solution 
only for > 1. To clarify the situation further, let us use equation (127) and obtain the 
dependence of r^^.^, ^ on the value of Tg/T. Substituting from equation (127) into equation 
(125) yields 
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Te//,. = 5.64/3^/^ (^^) (l-Co)-' (128) 

As /3 << 1, for Tg/T = .75 the optical depth exceeds unity only if i^o is close to one, which 
in turn by equation (127) requires > 1. Alternatively, substituting (1 — i^o) from equation 
(127) into equation (125) results in 



r.f,. = i.6fl"'^ (I)""" i.f/-Af''-(i -,»)-»'■■ (^y"'\iy^-^- (^)""". 

(129) 

leading to the same conclusion that > 1 is required. Moreover, substituting (1 — i^o) from 
equation (127) into equation (51) yields 



, /r< \ 8/33 /rp \ 20/11 

I = 0.34i:y- (1 - /3)^/- M-^/-eo-"/"r^/^^ (f^) (I) (130) 

which for r = 20 and Tg/T = 0.75 is ~ 0.7, so the disk is no longer geometrically thin. We 
conclude that in order for the modified blackbody approximation to be selfconsistent in this 
region, Tg/T must be smaller. In the absence of a detailed vertical structure, the actual value 
is not known. To illustrate the effect of a smaller Tg/T we adopt a value smaller by a factor 
of 10: 0.075. In this case equation (127) becomes 

(1 - eo)eo"/" = 0.03M-^/-i:;-/-(l - /3)-/n (^^^ ^l_y'^'' (131) 
and equations (129) and (130) yield 



r.„,, = 660fl'/'^if/ »A4"" (1 - ,»)-""' (1^) d"'" (^)""" ■ (132) 

A = omlV" (1 - 6f'" M-""/;;"''''' (^)°'" (|it)"" (133) 

For r = 20 equation (131) has a solution if > 0.1. The value of i^o ranges from ~ 0.45 
up to ~ 0.97 when varies from ~ 0.1 to ~ 1, respectively. For all these solutions the 
surface layer is optically thick for absorption (modified by electron scattering) and the disk 
is geometrically thin, even for > 1. 
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Substituting (1 — i^o) from equation (131) into equations (52)-(55) results in 



, = (3.5 X 10-^ gem--) M;'^ ' lI^ ' (1 - (^) (^) , (134) 

re. = 4.6 X lO^Mj/^L-/- (1 - /3)-/- il"' (^) (^) , (135) 

T = (2.9 X 10^ JO M-^Z-L^/- (1 - /3)^/- ^0"'/" (^) , (136) 

^ = 12Mj/-i:-/- (1 - /3)-/- Co""/" (^)''" (^)''" , (137) 

Note that the r dependences in equations (134)-(137) are quite different from those in equa- 
tions (52)-(55). The differences resuh from the substitution of (1 — 1^0) from equation (131) 
into the latter equations. In particular, with 1^0 given by equation (131), the ratio Pr/Pg 
increases with r since the dependence on r through 1^0 in equation (55) overcomes the explicit 
r dependence in that equation. Equation (137) implies that Pr > Pg for r > 6, provided that 
equation (131) is satisfied. Equation (131) has a solution for 

Thus for any given there is a maximal radial extent of the radiation pressure region 

/ 7- \ 64/93 X -16/93 

r<20(^) Afr(^) . (139) 

From equation (131) and (135) follows that (at fixed r) an increase of causes an 
increase in Tg.. Therefore, the radiation pressure region is secularily stable. The physical 
explanation for this stability is rather simple. Equation (53) indicates that Tg. is inversely 
proportional to and to (1 — i^o)^- From equation (131) follows that (1 — 1^0) -L*^^'^^^, 
thus the resulting increase in Tg. overcomes the decrease due to increase of L^. 
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The line marked R in figure 1 shows the luminosity versus the electron scattering optical 
depth of the disk (proportional to the disk surface density), in the gas pressure region for 
r = 20, /3 = 0.2, Ms = 1 and Chf = 30. The plot represents equation (135) with i^o obtained 
from equation (131). The positive slope indicates that the solution is secularily stable, in 
contrast with the prediction of the a disk model. 

Discussion 

This paper presents a self-consistent solution for thin accretion disks with turbulent 
convection. The energy release is due to the interaction of turbulent convection with the Ke- 
plerian shear via a turbulent viscosity. Since turbulent convection provides the disk viscosity, 
the convective flux (in addition to the radiative flux) could be important in the transport of 
the generated energy to the disk surface. Employing a model for turbulence, the turbulent 
viscosity and the convective flux were obtained as functions of the physical parameters of the 
disk, which in turn are controlled by the former two. Having a model for turbulence, there is 
no need to resort to phenomenological parameterizations of either the viscosity (as done in 
the a-disk models) or the convective flux (as done in the mixing length approach). 

Solutions for both, radiation pressure dominated (inner) and for gas pressure dominated 
(outer) regions, with either electron scattering or free-free absorption opacities, were obtained. 
They provide the midplane temperature, density, disk half-width and the ratio between the 
convective and total fluxes at a given distance, for given central mass and given accretion 
rate. The solutions in the various regions of the disk depend on i^o — the fraction of the 
convective flux out of the total flux. We use an approximate estimate for the ratio of the 
surface to midplane temperatures as well as expressions for the flux emerging from the disk 
top , in terms of T^. This yields an estimate for i^o as function of the local physical conditions, 
in the different regions. 

In the gas pressure dominated regions we flnd two distinct types of solutions. The flrst 
(corresponding to strong convection) is characterized by a convective flux much larger than 
the radiative flux. In this case, the blackbody region encompasses the entire gas pressure 
region and could as well extend down to the inner disk boundary, with no radiation pressure 
region. In the second class of solutions (corresponding to moderate convection), the convective 
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flux is of the order of the radiative flux (~ 1/3 of total flux) and aU three disk-regions would 
exist. In the radiation pressure dominated region i^o increases with and decreases with r. 
Thus, at r = 20, it ranges from ~ 0.5 when ~ 0.1 up to very close to 1 when > 1. The 
radial extent of the radiation pressure region is larger the larger is the accretion rate. The 
radiation pressure region is optically thick and geometrically thin even for super Eddington 
luminosities. 

All the disk solutions, including the radiation pressure dominated region, are secularily 
stable. In the radiation pressure region, an increase in the accretion rate causes an increase 
in the fraction of the convective flux. As a result the net dependence of the surface density 
on the accretion rate is through a positive power. The present solutions differ markedly from 
the a disk behavior, also in additional aspects. In the case of a disks, the radiation pressure 
dominated region becomes optically thin and geometrically thick for luminosities close to 
the Eddington luminosity. In the present case, all the solutions are geometrically thin and 
optically thick, even for luminosities equal to or exceeding the Eddington luminosity (see Fig. 
1). Therefore, for luminosities of the order of the Eddington luminosity, the local spectrum 
will be a modifled blackbody while in a disks it will be a bremsstrhalung or compotonized 
bremsstrhalung. 

It is interesting to note that the trend of convection to result in larger surface density, 
lower temperature and increased stability, is evident even in models where the disk viscosity is 
described by an a parameter and the convective flux by the mixing length approach (Milsom 
et al. 1994) 

We found that, for the same central mass and accretion rate, the innermost regions 
could be either radiation pressure dominated or gas pressure dominated blackbody with the 
convective flux much larger than the radiative flux. While the values of the corresponding 
effective a-parameter differ only by a factor of ~ 2 (see below), in the former case the disk 
is much thicker geometrically, much less dense and hotter than in the second case. The 
emerging flux is a modifled blackbody while in the second case it is a blackbody. Since most 
of the disk emission comes from the innermost regions, the spectral signature of these two 
solutions will be quite different for the same value of bolometric luminosity, corresponding to 
the same accretion rate. 
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As discussed above, the present disk model differs markedly from a disk models in various 
aspects. Nevertheless, it seems that the direct observational distinction between the two is 
not a straightforward one. For a given accretion rate, the temperatures and densities in the 
optically thick regions of an a disk are of a similar order of magnitude as in the present 
model, thus leading to similar spectra. As mentioned above, for high luminosities the spectra 
in the radiation pressure region would be quite different: a bremssrahlung or compotonized 
bremsstrhalung spectrum in the case of alpha disk versus a modified blackbody in the present 
case. However, an optically thin corona on top of an optically thick disk would also produce 
a hot bremsstrahlung or compotonized bremsstrhalung spectrum (Haardt & Maraschi 1991, 
1993). 

The values of the effective a-parameter in the different regions are quite small: < 5 X 
10""^, ~ 1 X 10~^, and ~ 5 X 10~^ for radiation pressure region and for the two types of 
solutions in the gas pressure region, respectively. These small values reflect the fact that the 
turbulent velocities are subsonic and that the horizontal scale of the largest eddies is small 
compared to the vertical scale. This anisotropy, resulting from the Coriolis force, is a general 
feature of a three-dimensional turbulence in a rotating disk. The above values are small 
compared to values required to model outburst of cataclysmic variables (see, e.g., Duschl 
1989). This implies that turbulent convection cannot account for the disk viscosity in this 
case. A different mechanism, possibly magnetic viscosity, is required. One may speculate 
that turbulent convection can be the source of the disk viscosity during the quiescent state 
and some other mechanism provides the disk viscosity during the outburst. 

Finally, we wish to stress the need for a model capable of providing a detailed vertical 
structure yet being self consistent. This is no easy task as present available spectral turbulence 
models, like the one applied here, assume homogeneity in the vertical direction and thus 
can yield only vertically averaged description. Also, the incorporation of the interaction of 
shear with the eddies is an open question in such models. A promising alternative is the 
Reynold stress approach that does not provide the turbulence spectral function but rather 
supplies differential equations in space and time for the various turbulence ensemble averages. 
The formalism has been employed, extensively and quite successfully, in atmospheric and 
laboratory turbulence (Zeman 1981; Speziale 1991) and has been used recently to study 
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convective overshooting (Canuto 1992, 1993) in stars. 
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Log ( Electron Scattering Optical Depth) 
Fig. 1. — Lminosity, in units of the Eddington luminosity vs. the electron-scattering optical 
depth, for the various convective disk solutions, at a constant distance r = 20, for 
Ms = 1 and Chf = 30. R: radiation pressure region, fi = 0.2, eqs. (135), (131). BBS: 
gas pressure blackbody solution in the strong convection limit, eq. (93). BBM: same 
region for the case of moderate convection, eq. (109). ESM: gas pressure electron- 
scattering solution in the moderate convection limit, eq. (117). 
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